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We study the interaction of Anderson localized states in an open ID random system by varying 
the internal structure of the sample. As the frequencies of two states come close, they are trans- 
formed into multiply-peaked quasi-extended modes. Level repulsion is observed experimentally and 
explained within a model of coupled resonators. The spectral and spatial evolution of the coupled 
modes is described in terms of the coupling coefficient and Q-factors of resonators. 
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Transport in open disordered media can be diffusive 
or localized, depending on the nature of the underly- 
ing quasimodes, which are, respectively, spread through- 
out the sample or exponentially peaked within the sam- 
ple, with a typical size given by the localization length 
HI, 0- Q • The spatial overlap of localized modes which are 
close in frequency, couples these states and leads to the 
formation of a series of exponential peaks known as neck- 
lace states [4, 5, 6|,|7|,l8|- These states are short-lived with 
broadened spectral lines 0, 0] and contribute substan- 
tially to the overall transmission in samples much thicker 
than the localization length, L ^> l\ oc [y, [8[. Though 
such hybridized states are critically important in trans- 
port and play a significant role in the localization transi- 
tion, their formation and the correlation between spatial 
and spectral properties has not been explored. 

In this Letter, we study the transformation of cou- 
pled Anderson localized states in a random sample as its 
configuration is altered leading to the hybridization of 
modes. Although level repulsion is ordinarily associated 
with the diffusive regime [10|, UjJ, energy level correla- 
tion and repulsion in localized ID electro n sy stems has 
been found theoretically and numerically 12, 13j]. Here 
we present the first experimental evidence of level repul- 
sion of localized electromagnetic excitations. A simple 
theoretical model is introduced which explains the spec- 
tral and spatial characteristics of modes in terms of losses 
within the sample and the strength of coupling between 
the modes. 

The experiment involves a rectangular microwave 
waveguide opened at both ends, which supports only a 
single transverse mode jgj. The waveguide is filled with 
a sample comprised of five 4 mm-thick blocks each of 
low and high indices of refraction randomly mixed with 
31 randomly oriented 8 mm-thick binary blocks with low 
and high index halves. The field inside the sample is 
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FIG. 1: (Color online.) Resonant frequencies of excited local- 
ized modes vs. the driving parameter - the air gap inside the 
sample. Five pair-interaction regions are circled. 



weakly coupled to a cable which is translated along a 2- 
mm-wide slot cut along the waveguide in 1 mm steps. 
Field spectra are measured using a vector network ana- 
lyzer. 

Measurements are made in a sequence of configura- 
tions in which the spacing between two segments of the 
sample at a depth of 60 mm is increased in steps of 0.5 
mm up to a maximum thickness of 14 mm. The position 
at which the air gap is introduced was chosen to corre- 
spond to the peak of a single Anderson localized mode 
of the unperturbed random sample. This allowed us to 
tune the frequency of the selected mode in a manner sim- 
ilar to the tuning of a defect states through a band gap 
in a periodic structure. However, here we deal with the 
Anderson localized states arising from the interference 
of the multiply scattered fields in a statistically homoge- 
neous random system. The mode frequency shifts across 
the frequencies of other localized states which makes it 
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FIG. 2: (Color online.) Experimentally measured normalized 
wave intensity vs. frequency and coordinate inside the sample 
for two localized modes corresponding to region 2 in Fig. 1 at 
different values of the driving parameter. Excitations corre- 
sponding to the two effective resonators are circled. 



possible to study the coupling of two modes. 

The spectral positions of the localized states as func- 
tions of the air gap introduced into the sample are plotted 
in Fig. 1. The frequencies of modes may cross or exhibit 
an anti-crossing (level repulsion). In the latter case (re- 
gions 1,2,4,5 of Fig. 1), the coupling within the sample 
is accompanied by an exchange of shape, as is seen in 
Fig. 2. When the frequencies of the modes are closest, 
the two localized states couple into double peaked quasi- 
extended modes with the same spatial intensity distribu- 
tion, Fig. 2b and c. In contrast, region 3 in Fig. 1 shows 
a mode crossing in which the shapes are not exchanged. 
This is seen in Fig. 3 which shows the driven mode pass- 
ing through the broad mode closest to the input. The 
two modes remain practically independent of each other, 
except for the low-intensity zone (dark horizontal line in 
Fig. 3) at the input mode. 

Resonant wave transmission through an isolated local- 
ized state in a random sample can be described by a 
simple model of a wave tunnelin g th rough a resonator 
with semitransparent walls 

BQHii- Dynamics of the 
field in the resonator obeys the oscillator equation with 
an external force and damping, which accounts for the 
incident wave and the finite Q-factor of the resonator, 
respectively. Extending this model to the case of N lo- 
calized states which are close in frequency, we arrive at 
a system of N coupled oscillators with the external force 



FIG. 3: (Color online.) The same as in Fig. 2 for interaction 
region 3 in Fig. 1. 



acting on the first of these: 



+ + (i - A;) 2 ipi = ©i+iVi+i + qu-ii>i- 



k 4>N + QjvVn + (1 - An) 2 ^n = qNN-li> 



JV-1 • 



Here ipi (r) is the field in the ith resonator, r = coot 
is the dimensionless time (coo is a characteristic central 
frequency of the problem), 1 — (|A^| <C 1) is the di- 
mensionless eigenfrequency of ith resonator, Qi ^> 1 is 
the Q-factor describing the losses of the energy in the ith 
resonator, qa+i = qi+n <C 1 is the coupling coefficient 
of ith and (i + l)th resonators due to the spatial overlap 
of their modes; /o and ^, (\v — 1| <C 1) are the ampli- 
tude and frequency of the external field exci ting the first 
resonator. The Q-factors can be written as [161 ]: 



r, (i< i < n) , 



'l,iV 



N ' 



, (2) 



where Ti is the dissipation rate in the ith resonator, 
Ti n , out are the transmission coefficient of the input and 
output of the system, v g is the wave group velocity inside 
the resonator cavity, and I is the cavity length. The last 
term in Eqs. (2) accounts for the energy leakage through 
the outermost resonators. 

To establish the correspondence between the model 
Eq. (1) and localized states in a random sample, we as- 
sume, following 2IIH, that ^ represents the peak field 
of the ith localized state, qa+iijji is the amplitude of the 
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FIG. 4: (Color online.) Field amplitudes |^4i,2| 2 in the two 
resonators as functions of the incident field frequency v and 
detuning A between the resonators. The underlying frequen- 
cies 5is^l s are depicted by the dashed lines. Parameters are: 
f = 1, q = 0.2, and Q" 1 = 0.1 (Q _1 < q regime). 



field penetrating the adjacent cavity, and fo is the am- 
plitude of the incident wave, ^o, penetrating the first 
localization cavity. Since for long enough system the lo- 
calization length is the only disorder-induced spatial scale 
in the problem, we assume 



q ii+1 ~ exp(-d ii+1 /li oc ) , f ~ ip exp (-d in /li oc ) 
— exp ( 2>di n ^ ou t I lioc) i I ^ hoc • 1 



i, out 



(3) 



Here da+i = — Xi\ is the distance between adja- 

cent states at coordinates Xi+\ and Xj, whereas di n = 
X\ and d out = L — Xn are the distances from the edge 
localized states to the corresponding ends of the sam- 
ple. The localization length is determined through the 
average value of the total transmission coefficient T as 
l~l = - < lnT(L) > /2L [3]. The deterministic equa- 
tions (l)-(3) provide an effective description of coupled 
modes in ID random system. 

Substituting ^ = Ai exp ( —ivr), the set of equations 
(1) is reduced to an algebraic equation = F with 
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^ = (A U ...,A N ) T , and F = f (1, 0, .., 0) T , where 
C % = (1 - A,) 2 - v 2 - ivQ~\ ^ 2 (1 - A, - v) - iQ-\ 
The homogeneous equation = determines a set of 
eigenmodes, with eigenfrequencies being the eigenvalues 
of the matrix (4). 

For the sake of simplicity, we consider the case of two 
interacting modes, N = 2, and assume that Qi = Q2 = 
Q and qi2 = q. Then, the complex eigenfrequencies z^ ± = 
1 + Sv ± are given by 

Ai + A 2 
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FIG. 5: (Color online.) The same as in Fig. 4 for the param- 
eters q = 0.2 and Q" 1 = 0.25 (Q _1 > q regime). 



This equation describes level anti-crossing and the cou- 
pling of isolated resonators to form collective eigenmodes. 
The minimal frequency gap q is achieved at resonance, 
Ai = A2. Away from resonance, |Ai — A2I ^> q, the 
eigenmodes tend to the modes of isolated resonators. The 
shapes of the modes are exchanged when passing through 
the resonance, i.e., + (— ) eigenmodes correspond to the 
first (second) resonator at Ai <C A2 and to the second 
(first) resonator when Ai ^> A2. It is important to note 
that level repulsion of electromagnetic modes arises in a 
finite system in the regime of strong localization. 

If the system is excited by a monochromatic wave with 
real frequency v = 1 + as in the experiment, the 
complex amplitudes Ai^ of the two resonators can be 
obtained from = F, which yields 
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[2 (Ai + Sv) + iQ^ 1 } [2 (A 2 + Sv) + iQ? 1 ] - q 2 

The behavior of |^4i, 2 | 2 is essentially determined by the 
denominator \D\ 2 , which is at minimum at frequencies 
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The amplitudes A\^ and frequencies Sis^ s characterize 
the resonant excitation of the system by an external 
source. Note that Eq. (7) coincide with Eq. (5) only 
in the lossless case Q _1 = 0. Otherwise, there are two 
different regimes of the excitation of coupled resonators, 
determined by the ratio between losses Q _1 and coupling 
q. If losses are small, Q~ x < q, two branches Sv^ s demon- 
strate anti-crossing with a frequency gap y / 'q 2 — Q~ 2 : 
Fig. 4. If losses are greater than the coupling strength, 
so that Q _1 > q, the frequencies 5u^ merge in the inter- 
action region (Ai - A 2 ) 2 < \q 2 - Q _2 |, Fig. 5. 

The amplitudes (6) for Q _1 < q and Q~ l > q are 
shown as functions of the frequency v and detuning 
(Ai — A 2 ) in Figs. 4 and 5, respectively. The main fea- 
tures observed experimentally are reproduced. Different 
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FIG. 6: (Color online, 
retically calculated (Ke^/q 2 — Q 
for pairs of interacting modes 1,2,4,5 presented in Fig. 1 



Experimentally measured and theo- 
2 ) minimal frequency gaps 



values of A correspond to different frames of Figs. 2 and 
3. To facilitate comparison with the experimental results 
in Figs. 2 and 3, the second, output resonator is driven 
in Fig. 4 (Ai =0, A2 = A), while the first resonator is 
driven in Fig. 5 (Ai = A, A2 = 0)). It is seen in Figs. 
2 and 4 that for Q~ x < q fields in both resonators in 
the interaction regime (Figs. 2b and c) exhibit double- 
peaked spectra (level repulsion). Collective excitation of 
two resonators signifies the formation of quasi-extended 
necklace states. Remarkably, away from the resonance 
(Figs. 2a and d), the first resonator is effectively excited 
at one of the resonant frequencies, close to Su = — Ai 
(Fig. 4a), while the second resonator is equally excited 
at both resonant frequencies Sv^ s ~ Ai } 2 (Fig. 4b). In 
the regime Q _1 > <?, both Figs. 3 and 5 show that the 
second resonator is excited with a single peak in the spec- 
trum (Fig. 5b), while two peaks separated by a dark area 
driven with the frequency of the second resonator are seen 
in the first resonator (Fig. 5a) [17j . 

The measured and calculated values of the frequency 
gap between coupled modes are presented in Fig. 6. The 
parameters of the system are: 



^ ~ 15.5GHz, ho 



12mm, cj r - 7x 10' s 



7.-1 



(8) 



and v g ~ c/2.4, whereas the positions of the localized 
modes interacting in the regions 1-5 (Fig. 1) are equal, 
respectively, to 



x 2 



64 : 64 : 7 : 64 : 64 mm , 
117 : 192 : 64 : 128 : 235 mm . (9) 



Substituting values (8) and (9) into Eqs. (2) and (3) 
yields Qi \ and q. We calculated the minimum frequency 
difference for interacting pairs 1,2,4,5 (for which Q^ 1 ~ 
Q2 1 [17]) as Re^q 2 - Q~ 2 with Q' 1 = (Qf 1 + Q 2 " 1 )/2, 
and compared it to measured values of the gap, Fig. 1. 
Good agreement between the experiment and model is 
seen in Fig. 6. 



In conclusion, we have observed level repulsion in the 
localization regime and have shown that it reflects the 
coupling of localization centers. The occurrence of anti- 
crossing or crossing of quasimodes as a sample configu- 
ration changes depends upon the ratio of the coupling 
strength between localized states and loss. These factors 
determine the statistics of level spacings and widths and 
thus the diffusive or localized regime of wave propaga- 
tion. 
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It should be noted that two localized modes of Fig. 3 
have significantly different Q-factors, Q^ 1 ^> Q^ 1 • This 
diminishes the coupling between the modes since their 
frequencies are distant in the complex plane. Nonetheless, 
our simplified model with Q\ — Q2 describes quantita- 
tively all the main features of the coupling in this case. 



